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Abstract. We have investigated the coherent phase control on the autoionizing state (AIS) resonantly 
coupled with the ground state for Mg through a two- and a four-photon transition simultaneously, using a 
bichromatic linearly polarized laser field. The frequency is chosen such that the lasers are tunable around 
resonance with the transition 3s^(^5"') 3p^(^5"'), which implies ojf — 2.11 eV and ujh = 4.22 eV. We 
are interested in the modification of autoionizing (AI) line shape through the relative phase and laser 
intensities. A strong phase dependence on the total ionization yield and ionization rate is found. We 
also performed a time-dependent calculation which takes into consideration all the resonant states of the 
process. 

PACS. 3 2.80.Qk Coherent control of atomic interactions with photons - 32.80.Rm Multiphoton ionization 
and excitation to highly excited states (e.g., Rydberg states) - 32.80.Dz Autoionization 



1 Introduction 

This paper treats a particular aspect of the control of pho- 
toabsorption, and in particular, ionization through the rel- 
ative phase of two electromagnetic fields acting simultane- 
ously on an atomic system. The aspect in question is the 
possibility of altering the line shape of an autoionizing res- 
onance through such phase control. The general features 
of phase control in photoabsorption have been discussed 
in the literature extensively [1,2, 3]|31[SJini[Z] • The possibil- 
ity of altering an autoionizing line shape has also been 
discussed [5] to some extent in theoretical work where it 
has been shown that one particularly intriguing feature, 
depending on parameters, is the possibility of turning off 
the transition to the "discrete" or the continuum part of 
the resonance [51ITU]. To the best of our knowledge, experi- 
mental observation of such an effect has not been recorded, 
one of the reasons perhaps being the unavailability of ra- 
diation sources of intensity and frequency suitable for con- 
venient atomic species and resonances. 

The case studied in this paper aims at presenting a 
quantitative analysis of a situation possibly convenient 
for existing laser sources whose potential in phase control 
has been tested experimentally. We have thus undertaken 
the study of the line shape of the 3p^ i^S'^) autoionizing 
resonance of atomic Magnesium under the simultaneous 
excitation by four- and two-photon transitions whose rel- 
ative phase is varied. This case [TT] is somewhat different 
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from, and theoretically a bit more demanding, than the 
more standard scheme of one- and three-photon transi- 
tions, which in any case would not be applicable here ow- 
ing to the parity of the resonance which is the same as that 
of the ground state. In addition, we have found that bound 
states in near resonance with one- or three-photon transi- 
tions introduce distortions in the wings of the autoionizing 
resonance; an effect not a prior obvious but which has to 
be taken into account for a realistic assessment of the ob- 
servability of the desired feature. Needles to say that a 
realistic atomic structure and transitions calculation has 
been necessary, whose details are discussed in the sections 
that follow. 

Lyras and Bachau [8; studied the phase control in two- 
and four-photon ionization of the Magnesium atom in a 
bichromatic field of frequencies w and Sw, in the vicinity 
of an autoionizing resonance ^D^ lying above the first two 
ionization thresholds. Moreover, they accomplished a full 
perturbative, time-independent phase control calculation 
in Mg by interfering three- and one-photon transitions to 
a single continuum channel ^P". They studied the depen- 
dence of the ionization rate as a function of the relative 
phase for different laser intensities, and noticed that the 
presence of an intermediate resonance may enhance the 
overall ionization signal. 

Kylstra et al. [12] have performed a non-perturbative ab- 
initio one- and two-color calculation of the multiphoton 
ionization of Magnesium, where the laser frequencies are 
chosen such that the initial state of the atom is resonantly 
coupled to the autoionizing 3p^, 3s3d, and 4s5s resonances 
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of the atom. Using R-matrix Floquet theory, they studied 
single photon ionization from the ground state 3s^ and 
the 3s3p excited state of Magnesium in the vicinity of the 
AI resonance at non-perturbative laser intensities. A 
second low-frequency laser couples the and 3s3d AI 
states. This resonant system was studied experimentally 
and theoretically in the perturbative domain. Kylstra et 
al. also investigated the harmonic intensity regime where 
the process is no longer perturbative. 

In section 2 we present the basic framework related 
to the study of the phase control in the vicinity of an AI 
state. In section 3 we briefly give the ab-initio theoreti- 
cal approach for the description of the Magnesium atomic 
structure. Further details can be found in P^fTHITS] . Fi- 
nally, in section 4 we present our results. 



2 Basic Equations 

We study the multiphoton ionization of Mg using a realis- 
tic atomic model which describes the time-evolution of the 
system, investigating the phase coherent effect on the 
i^S'^) autoionizing state resonantly coupled to the ground 
state of Mg, simultaneously through two- and four-photon 
transitions, in the weak field regime. The laser frequency is 
such that the ground state of the atom is near resonance 
with the 3s3p {^P") {EssSp ~ 4.34 eV) state of the Mg 
atom due to its second harmonic, and with 3s4p C'P") 
{Essip ~ 6.11 eV) due to three-photon absorption from 
the fundamental. 

We consider the Magnesium atom consisting of a ground 
state \g >= |3s^ ^S"^ >, an autoionizing state \a >= |3p^ 
^S'^ >, and one continuum corresponding to three dif- 
ferent values of angular momentum |ci >= |3ses ^S'^ >, 
|c2 >= |3sed iD'^ >, and |c3 >= |3se<? ^G^ >. The AIS is 
modeled as a discrete state embedded into the continuum 
and coupled to the ^ S'^ continuum through the configura- 
tion interaction. 

In Figure [1] we show the energy diagram of the atomic 
system interacting with the linearly polarized bichromatic 
field: 

E(i) ~ E/(t) exp {iujft) + 'Eh{t) exp {iujht + iip) + cc., (1) 

consisting of a superposition between the fundamental 
with frequency uj f and its second harmonic with frequency 
Ldh = 2a;/. The amplitudes of the two components of the 
electric field are £f{t), £h{t), and (p is the relative phase 
between them. The continuum state of energy Eg + Aujf 
can be reached by the two interfering paths shown in Fig- 
ure [1] namely, the four photon absorption from the funda- 
mental electromagnetic field (Figure [1] path (i)), and the 
two-photon absorption of its second harmonic (Figure [1] 
path (ii)). 

There are a few other processes which might affect the 
interference process, but for the intensities we consider in 
this paper they are not significant, therefore we do not 
need to include them in the calculation. For instance, a 
three-photon transition with two photons of the funda- 
mental and one of the harmonic could lead to the same 



final continuum energy. These processes, however, could 
only influence the background signal, and not the interfer- 
ence process, since the final continuum state belongs either 
to a ^P" or to a state. We have estimated the transi- 
tion amplitudes for these three-photon transitions and for 
values of the laser intensities involved in our calculation, 
and we found out that the contribution to the ionization 
signal is about two orders of magnitude smaller than the 
contribution due to the processes presented in Figure [1] 
One can therefore conclude that these three-photon pro- 
cesses do not affect the process studied here and we neglect 
them. The dominant processes are decided by the value of 
the frequency and intensities involved in the calculation. 
For the same reason, Raman-type processes are ignored 
due to the fact that they represent higher-order processes 
with respect to the electric field. We have used atomic 
units throughout this work. The atomic unit used for the 
intensity of laser field is /o=14.037 x 10^^ W/cm^. 

A four photon transition from the |3s^ ^S"^ > state into 
the continuum leads to a final state containing \ 3ses ^S"^ >, 
\3sed ^D'^ >, and \3seg ^G^ > partial waves, while a two- 
photon transition into the continuum leads only to |3ses ^S*"^ > 
and |3se(i > partial waves. Thus, only the transi- 
tion amplitude to the ^S'^ and ^D'^ continua is modu- 
lated through the interference with the two-photon ampli- 
tude. The transition amplitude to the continuum re- 
mains unaffected by the quantum interference, and there- 
fore only contributes to the background of the total ion- 
ization signal. 

We consider the Schrodinger equation: 

i^^[H,+D{m{t), (2) 

where Ha represents the atomic Hamiltonian. The atomic 
system is treated as a two active electron system: the 
atomic core (the nucleus and the 10 inner-shell electrons), 
and two valence electrons. More details about the atomic 
structure calculation are given in the next section. Within 
the semiclassical formalism, the interaction between the 
atom and the bichromatic field is described in the length 
gauge, and in the dipole approximation by D{t) — — (ri -|- 
Y2) ■ E(i), where ri and Y2 represent the position vectors 
of the two valence electrons. In order to describe the dy- 
namics of the system, the time dependent wavefunction 
is expanded in terms of the complete set of states, and 
then substituted into the Schrodinger equation following 
the standard procedure [TSJ. Pd t — {) the system is as- 
sumed to be in the ground state \g >, and for time i > 
its wave function can be written as: 

3 

where Cg, and Cj are the probability amplitudes of states 
\g > and \cj >, assumed to be eigenstates of the atomic 
Hamiltonian Ha with the eigenenergies Eg and Ec- , re- 
spectively. The substitution of the wave function Eq. 
into the Schrodinger Eq. ([2|) leads to a set of differen- 
tial equations for the amplitude coefficients. We intro- 
duce the slowly varying coefficients Cg = CgC*^^"-'* and 
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Cj = Cje''-'^f+^"^'*, and, after the adiabatic elimination of 
the coefficients for the continuum \cj > in the differential 
equation for the coefficient of the ground state, using the 
rotating wave approximation, we are left with one inde- 
pendent equation (for the ground state) , and three coupled 
differential equations governing the time evolution of the 
probability amplitudes as such: 



. dCgjt) _ 
' dt " 

,dcj{t,Ec. 



dt 



(4) 
(5) 



+ 



.dc3{t,Ecs) 



dt 



C3{t,E,,)e,,+DifJE,„t)cg{t), (6) 



Er_ 



where - j^c^ - ^g 
We used the identity: 

hm 

17^0+ J X — xo + ir] 



Eg-4ojf, i = l, 2, 3. 



dx^V 



X — Xo 



-dx - infixo), (7) 



where V is the Cauchy principal value integral. Sg repre- 
sents the ground state AC- Stark shift induced by the laser 
field. Since we are working in the weak laser field regime, 
Sg wiU be neglected. 

The quantities Dcf'g and oi^^l are the effective two- and 
four-photon transition amplitudes from the ground state 
\g > to continuum \cj >: 



n 



J = 1,2 



(8) 



Mif2 Af/f ikfi^i 



J = 1,2,3, 



Jng - iu:f){uilg - 2U} f){u},ng - ^ f) 

(9) 



where the one-photon transition amplitude between the 
state I A > and \B > are written in the length gauge as 
^abW = - < A|(ri+r2)-E^(0|S >, andMy](i) = - < 
A|(ri -|-r2) •E/(i)|i? >, respectively, and ujab = Ea — Eb- 
The quantity 7^ represents the effective total ionization 
width directly into the continuum, from the ground state: 



I9 



2iip 



-2tt 



Dm 



(10) 

where the bar over the effective transition amplitudes Dgc - , 

and D^'^J. means that they have been calculated at Ec — 
Eg + iiOf,j ^1,2,3. 

Firstly, we consider the ionization probability per unit 
time P, which is valid for an almost square pulse £f = 
£f{t), and £h = £h{t)- In this case, the photoionization 
line shape is simply obtained in terms of the transition 
rate without any time-dependent calculations by assum- 
ing dtCj{t) = 0. The ionization rate is given as the loss of 



population from the ground state: 

P = ~2lm[zdtCg{t)c;{t)]^jg{^). 



(11) 



It is possible to control different ionization products through 
the relative phase of the laser components. The branching 
ratio into the j channel is defined as Bj = Pj/P, where 
the partial ionization rates Pj are: 



P:, 



27r 
27r 



1,2 



(12) 
(13) 



Our calculation takes into account all the important 
resonances with intermediate atomic states for the pro- 
cess studied. For the range of intensities where the two- 
and four-photon transition amplitudes are of comparable 
values, quantum interference effects of the two ionization 
channels would determine the modulation of the ioniza- 
tion signal. The photoionization signal and the AI 
line profile can be controlled by varying the relative phase 
ip between the field components of the laser field E(t). 
If we need to explicitly take into account the temporal 
evolution of the laser pulse, we have to integrate the time- 
dependent system of differential equations Eqs. (jH)-®, 
with the initial conditions Cg(i = 0) = 1 and Cj{t = 0) = 
0, j = 1, 2, 3. The temporal pulse shape of the laser field 
is a sine-squared function, £i{t) — £i sin^ {nt /T) , where 
i = f,h. The integration time for this sine squared shape 
pulse is taken from t = Q to t = T 

the number of laser cycles. 

Additional parameters, such as the laser pulse dura- 
tion and the temporal overlap of the fundamental and 
harmonic pulse, may also control the ionization signal. 

The total ionization probability is calculated at the 
end of the laser pulse: Pion{T) = 1 — |cg(T)p. We have 



— n, with n being 



obtained the value P- 



3p2 



1.76 X 10 a.u. for the au- 



toionization width of the 3p^ AIS, and the corresponding 
life-time of this level is about T^p2 = 13.75 fs. The total 
laser pulse duration used in our calculation is larger than 
the life-time of the AIS. 



3 Atomic structure 

The Hamiltonian of the Magnesium atomic system Ha is 
given by: 



Ha = 



2 

E 

1=1 



1 



(14) 

where V^^^(ri) is the radial Hartree-Fock potential for 
the closed-shell core of Magnesium (Mg"*""*"), and Vd is a 
two-body interaction operator that includes a dielectronic 
effective interaction P^ITTj . The characteristic feature of 
the Magnesium atom is the existence of a 3s^ valence shell, 
outside a closed-shell core, the excitation of which requires 
a much larger amount of energy compared with the first 
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and second ionization threshold. This allows us to explore 
excitation and/or ionization processes of the valence elec- 
trons, for a certain range of photon energies, without con- 
sidering the closed-shell core-excitation. 

At a first stage, we perform a Hartree-Fock calcula- 
tion for the closed-shell core of Magnesium (Mg++), thus 
deriving the effective Hartree-Fock potentials 'seen' by the 
outer electrons. At a second stage, we solve the Schrodinger 
equation for the valence electrons through the configura- 
tion interaction (CI) method using the Magnesium Hamil- 
tonian defined in Eq. ((T?)) . In this case, we also add a core 
polarization potential acting on the valence electrons, the 
inclusion of which represents the influence of the core on 
the two valence electrons in a way very similar to that de- 
scribed in [T^fTl] . This core-polarization potential has the 
form as {l — exp[— (r/r;)^]} /r'^, where as is the static po- 
larizability of the doubly-ionized Mg, and the cut-ofi^ radii, 
ri, for the various partial waves Z = 0, 1, 2, ... |17| . 

Having produced the Mg"*" one-electron radial eigen- 
states Pni{r) for each partial wave I = 0,1, 2, we solve 
the two-electron Schrodinger equation: 



(15) 



by expanding the two-electron eigenstates iJ'"^(ri,r2) on 
the basis of LS-uncoupled two-electron antisymmetrized 
orbitals ^;^;„,;,(ri, r2), namely 14 : 



C^in'i'{E,)<Pi^,i,{ri,r2), (16) 



iln'l' 



Pni{ri)Pn-v{r2) 



r2 



Xfi,f2;l,l'), 



(17) 

where A — (LM^), and A12 is the antisymmetrization 
operator which ensures that the total wave function is an- 
tisymmetric with respect to interchange of the space coor- 
dinates of the two electrons. Assuming the Magnesium is 
initially in its ground state |3s^ ^S"^ >, and that the tran- 
sitions with a linearly polarized light are well described in 
the dipole approximation, we only need to construct the 
singlet states with S = and Ml = 0. 

We force the wavefunction to be zero at the bound- 
aries by selecting the basis functions Pni(r), and Pn'i'{r) 
to be the one-electron radial solutions of Mg+ which, by 
construction vanish at the boundaries. The radial func- 
tions Pni(r) are expanded in a B-spline basis of order n, 
Pniir) = J2i biBi{r), i = 1,2, N, which transforms the 
one-electron Schrodinger equation for the Mg"*" orbitals 
into a system of matrix equations for the coefficients bi 

Substitution of the two electron wavefunction Eq. (jl7p 
into the Schrodinger equation Eq. ([T5]), leads to a gen- 
eralized eigenvalue matrix equation, the diagonalization 
of which gives the coefficients Cnin'i'{Ei) for each dis- 
crete eigenvalue Ei [T^lfTS] . This choice of the basis func- 
tions (ri, r2) (constructed from the radial orbitals 
Pni(r)), thus leads to a discretized continuum spectrum 
for the Magnesium atom with a density of states basically 
determined from the box radius R. 



Having obtained the two-electron wavefunctions '^nin'u ' 
we calculate the effective two- and four-photon transition 
amplitudes Eqs. (jS])-® within the lowest order-perturbation 
theory (LOPT), as well as the corresponding partial gen- 
eralized N- photon cross sections leading to the \cj > con- 
tinuum: 



aW(..) = 2.(2.a)^c.^ A^W( 



(18) 



where w represents the laser frequency, a the fine struc- 
ture constant, and j = 1, 2 for = 2 (two-photon cross 
section) and j = 1,2,3 for A^ = 4 (four-photon cross sec- 
tion). The effective two- and four-photons transition am- 
plitudes are calculated, D^g^^{uj) = I^^^Mg^^ (uj), where 
/ represents the laser field intensity. We have employed a 
box with the radius close to i? « 300 a.u., 602 B-spline 
functions of order nine, and a total number of angular 
momenta up to L = 4. 



4 Results and discussion 

In the first part of this section we present a few results 
in the weak field limit for a time- independent laser pulse, 
in which case the ionization rate approximation is valid, 
and the photoionization signal is well described by the 
ionization rate formula Eq. (jlip . 

Figure [2] shows the partial two-photon ionization cross 
sections from the ground state of the Magnesium atom 
leading to the ^S'^ (full curve) and ^D'^ continuum (dot- 
ted curve), as a function of the photon energy. Our results 
presented in Figure [5] are in agreement with the results 
published in the literature by Chang and Tang [19] and 
Kylstra et al. pjj. The four- photon partial cross sections 
from the ground state leading to the ^S^ (full curve), ^D'^ 
(dotted curve), and (dashed curve) continua are pre- 
sented in Figure [31 The two- and four-photon partial cross 
sections were calculated using Eq. 

The total ionization rate for the simultaneous two- and 
four-photon ionization from the ground state of Mg is plot- 
ted as a function of the fundamental laser frequency in 
Figure 01 for a relative phase 1^9 = 0° (dashed line) , 30° 
(dotted line), and 90° (full line). The laser field intensities 
are // = 2 x 10" W/cm^ and h = 4.13 x 10^ W/cm^. 
If the two fields are in phase, there is destructive interfer- 
ence between the two different ionization channels illus- 
trated in Figure [1] while there is constructive interference 
at (/3 = 90°. This behavior of the ionization rate is due 
to the relative sign of the effective two- and four-photon 
transition amplitudes; specifically, for a fundamental laser 
frequency at 2.08 eV< luj < 2.15 eV, the effective two- 
and four-photon transition amplitudes from the ground 
state into the ^S'^ and ^D"^ continua, respectively, have 
opposite signs. When the laser frequency is tuned near 
the resonant state 3s3p {ujh ~ 4.3, eV) no phase effects 
are observed, since the ionization through two harmonic 
photon absorption becomes dominant over the four pho- 
ton absorption. Figure [TJii), and no interference process 
exists. Where the laser is tuned around the 3s4p state 
{ujf w 2.02 eV), channel (i) of Figure [His dominant. 
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In order to have a better view of the interference pro- 
cess between the two ionization channels illustrated in Fig- 
ure m in Figure [S] we present the partial ionization rate 
into the ^S'^ (dot-dashed curve), ^D'^ (dotted curve) and 
^G^ continua (dashed curve). The full curve represents the 
total ionization rate when the fields are in phase, panel 

(a) , and when the relative phase is 1^9 = 90°, panel (b). 
We observe that the AI line profile for the two- and 
four- photon process is not symmetric due to the one- and 
three-photon transitions to the 3s3p and 3s4p near reso- 
nance states. 

Figure inja) shows the branching ratios function of the 
fundamental laser frequency aA, ip — 0°, and in panel b) 
for (fi = 90° with the same conditions as in Figure HI The 
branching ratio into the ^S"^ continuum is represented by 
the dot-dashed curve, into the ^D'^ continuum by the dot- 
ted curve and into the ^G"^ continuum by the dashed curve. 
As shown in panel (a), at (p = 0*^ for laser frequencies 
tuned around AIS, the ionization rate into the ^G^ 
continuum is enhanced since the laser intensities where 
chosen such that the two- and four-photon transition am- 
plitudes from the ground state into the ^S'^ continuum 
cancel each other. On the other hand, a,t (p — 90" in panel 

(b) there is an enhancement of the ionization signal into 
the ^S'^ and ^D'^ continua. This suggests that by a judi- 
cious choice of the laser intensities, it might be possible to 
arrange the maximum of one ionization product to coin- 
cide with the minimum of the others, and thus to control 
different ionization products. 

The modulation of the ionization rate as a function of 
the relative phase for a simultaneous two- and four-photon 
ionization from the ground state of Mg to the 3p^(^5^) 
AIS is depicted in Figure [71 for ujf = 2.11 eV with the 
same laser intensities as in Figure |4l The full curve repre- 
sents the total ionization rate, the dot-dashed curve is the 
contribution to ionization into the ^S"^ continuum, and the 
dotted curve the contribution coming from the ^D'^ contin- 
uum. The ionization rate into the continuum (dashed 
curve) is a flat line and contributes to the background 
of the ionization rate, which it cannot be neglected when 
the relative phase between the laser field components is 
close to kir. The ionization rate into the ^S"^ dX Lp = kir is 
zero since the laser intensities were chosen such that the 
effective transition amplitudes for two- and four-photon 
absorption into the continuum at w/ = 2.11 eV are 
almost equal, and cancel each other. 

We can analyze the ionization rate as a function of 
the harmonic intensity for fixed frequency and intensity 
of the fundamental. In Figure [5] we show the depth of 

ionization {Pmax - Prmn) / ^{Pmax + Pmin) aS a fuuctiou 

of the harmonic intensity. Pmax represents the maximal 
value of the ionization rate when <p = 90°, and Pmin is 
the minimal value of the ionization rate when the two 
laser components are in phase ip = O'^; the fundamental 
photon energy is w/ = 2.11 eV. As can be seen, an efficient 
coherent control in the weak field limit is obtained at If — 
2 X 10^^ W/cm^ for the harmonic intensity in the interval 
4 G (2 X lO^W/cm^, 2 X 10® W/cm^). Of course, good 



coherent control can be obtained as well for other pairs of 
laser intensities. 

When considering a temporal dependence of the laser 
pulse shape, the total ionization probability is obtained 
by numerically integrating the system of differential equa- 
tions for the amplitude coefficients in Eqs. (H))-®. In Fig- 
ure[9]we plot the ionization yield for field intensities where 
the depth of modulation has a maximal value: If = 2 x 
10^^ W/cm^ andJ/i = SxlO"" W/cm^. The ionization yield 
for the relative phase (/s = 0° is described by the dashed 
curve, for cp = 30° by the dotted curve, and for ip = 90° 
by the full curve. The laser pulse shape is sinusoidal with 
a total duration of 1 ps (n=500). 

Figures |3][7] and [5] clearly show that the AI line shape 
of is significantly changed by the relative phase ip 
of the laser components. The corresponding peak at 
ujf = 2.11 eV is diminished the most whenever the two 
fields are in or out of phase due to the destructive quan- 
tum interference between the two different channels, and 
enhanced the most whenever cp = {2k + l)7r/2 due to the 
constructive interference of the two ionization pathways. 
The maximum and the minimum of the ionization signal 
differ by more than one order of magnitude. 

In order to better illustrate the destructive interference 
between the two- and four-photon ionization processes, in 
Figure mH we plot the ionization yield as a function of the 
fundamental laser frequency at = 0°, and the same laser 
intensities as in Figure ^ The laser pulse shape is sinu- 
soidal with a total duration of about 1 ps. The dashed line 
represents the ionization rate corresponding to the four- 
photon absorption from the fundamental, and the dotted 
line describes the ionization yield due to the two-photon 
absorption of its second harmonic. The full line is the total 
ionization yield resulting from the destructive interference 
of these two processes. 

5 Conclusion 

In this paper, we have investigated the phase coherent 
effect of the 3p^ C'S'^) autoionizing state resonantly cou- 
pled to the ground state of a Mg atom in the presence of 
a bichromatic laser field of frequencies ujf and 2ujf. The 
transition amplitudes have been evaluated using a realis- 
tic atomic structure calculation. The motivation for this 
study was the investigation of the possibility of achieving 
coherent control of the photoelectron current and changes 
in the autoionization line shape. We have calculated and 
presented, for the first time (to the best of our knowl- 
edge), four-photon partial transition amplitudes for the 
Mg atomic system, and the corresponding partial ioniza- 
tion cross sections. We have found out that the near res- 
onance, double excited states 3s3p (^P") and 3s4p (^P°) 
of Magnesium introduce distortions in the left and in the 
right side of the autoionizing resonance, and enhance the 
overall ionization signal. A reliable study of the 3p^ line 
shape through the density matrix or resolvent operator 
method should take into consideration these two interme- 
diate resonances, and the energy dependence of the respec- 
tive transition amplitudes and Rabi frequencies. The rel- 
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ative phase ip between the two components of the bichro- 
matic laser field modulates the quantum interference be- 
tween the two ionization channels: four-photon absorption 
from the fundamental laser field and two-photon absorp- 
tion of its second harmonic. 

The work by G.B-Z. was supported by the European Research 
Network Program Contract No HPRN-CT-1999-00129. One of 
the authors (G.B-Z.) is indebted to Professor P. Lambropoulos 
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TABLE - FIGURE CAPTIONS 

Fig. 1. Energy diagram of the Mg atom interacting 
with a bichromatic electromagnetic field, showing the rel- 
evant levels for the studied process. 

Fig. 2. Two-photon ionization cross section cr'^-', lead- 
ing to (full line) and ^D^ (dotted line) continua from 
the ground state of Mg, as a function of the photoelectron 
energy. 

Fig. 3. Four-photon generalized ionization cross section 
cr(4\ leading to ^S" (full hue), ^D<' (dotted hne) and ^G" 
(dashed line) continua from the ground state of Mg, as a 
function of the photoelectron energy. 

Fig. 4. The ionization rate of the ground state of Mg as 
a function of the fundamental laser frequency for ip = 
(dashed line), 30° (dotted hne), and 90° (fuU line). The 
peak occurring at 2.02 eV corresponds to the 3s4p level, 
and the peak at 2.15 eV corresponds to the 3s3p level. The 
fundamental and harmonic laser intensity are If — 2x 10^^ 
W/cm^ and Ih = 4.13 x 10'' W/cm^, respectively. 

Fig. 5. The total and the partial ionization rates of the 
ground state of Mg as a function of the fundamental laser 
frequency. In panel (a) the relative phase between field 
components is 0°, and is 90° in panel (b). The partial ion- 
ization rate into the ^S*^ continuum is represented by the 
dot-dashed curve, into the ^D'^ continuum by the dotted 
curve and into the ^G"^ continuum by the dashed curve. 
The full line represents the total ionization rate. The laser 
parameters are the same as in Figure 21 

Fig. 6. Branching ratios as a function of the fundamen- 
tal laser frequency at (p = 0" in panel (a), and a,t (p — 90" 
in panel (b). The branching ratio into the ^S^ continuum 
is represented by the dot-dashed curve, into the ^D^ con- 
tinuum by the dotted curve, and into the continuum 
by the dashed curve. The laser parameters are the same 
as in Figure m 

Fig. 7. Ionization rate vs the relative phase between 
laser components. The laser intensities are the same as in 
Figure m The partial ionization rate into the ^S^ contin- 
uum is represented by the dot-dashed curve, into the ^D"^ 
continuum by the dotted curve and into the ^G"^ contin- 
uum by the dashed curve. The full line represents the total 
ionization rate. 

Fig. 8. The depth of modulation as a function of the 
harmonic intensity. The fundamental laser intensity is // = 
2 X 10" W/cm2. 

Fig. 9. The ionization yield of the ground state of Mg 
around the 3p^ AIS function of the fundamental laser fre- 
quency for = 0° (dashed line), 30° (dotted line), and 
90° (full line). The fundamental and harmonic peak laser 
intensity are // = 2 x 10^^ W/cm^ and Ih = 5 x 10'' 
W/cm^, respectively. The laser pulse shape is sinusoidal 
with a total duration of 1 ps. 

Fig. 10. The ionization yield at = 0° as a function of 
the fundamental laser frequency (full curve). The dotted 
curve represents the ionization yield for the two-photon 
absorption, and the dashed one corresponds to the four- 
photon absorption. The peak laser intensities are the same 
as in Figure [H 
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Fig. 10. 



